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ABSTRACT 

Following our previous work of Lachowicz & Czerny (2005), we explore further the 
application of the continuous wavelet transform to X-ray astronomical signals. Using 
the public archive of the XMM-Newton satellite, we analyze all available EPIC-pn 
observations for nearby Seyfert 1 galaxies MCG-6-30-15 and NGC 4051. We confine 
our analysis to 0.002-0.007 Hz frequency band in which, on the way of theoretically 
motivated premises, some quasi-periodic oscillations (QPOs) are expected to be found. 
Wc find that indeed wavelet power histogram analysis reveals such QPOs centered at 
two frequencies of ^2.5 x 10~ 3 Hz and ^4 6 x 10 -3 Hz, respectively. We show that 
these quasi-periodic features can be disentangled from the Poisson noise contamination 
level what is hardly to achieve with the standard Fourier analysis. Interestingly, we 
find some of them to be in 2:1 or 3:2 ratio. If real, our finding may be considered as a 
link between QPOs observed in AGN and kHz QPOs seen in X-ray binary systems. 
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1 INTRODUCTION 

One of the most remarkable features of the X-ray emission 
coming from active galactic nuclei (AGN) is their time vari- 
ability. Rapid and aperiodic, rather than predictable and de- 
terministic (e.g. Lawrence et al. 1987; McHardy & Czerny 
1987; Mushotzky et al. 1993), was quickly linked to processes 
taking place in the vicinity of the central black-holes. In a 
recent widely accepted unified model of AGN (e.g. Urry & 
Padovani 1995), soft X-rays/UV photons are emitted close 
to the central object in accretion disk and undergo upscat- 
tering via the inverse Compton process to higher energies 
in a hot plasma (e.g. Sunyaev & Titarchuk 1980). Thus, in- 
formation about intrinsic physical properties of the X-ray 
emitter is transmitted to the observed X-ray flux which is 
affected by the Poisson noise process occurring as a result 
of their detection (e.g. van der Klis 1989). 

The most classical way of description of broad-band 
variability of AGN is an estimation of the Fourier power 
spectral density (PSD). It gives an idea about the power of 
variable components (in terms of mean of the squared am- 
plitude) as a function of frequency or, inversely speaking, at 
characteristic time-scales related to the source. A commonly 
applied practice for all AGN is to study the overall shape of 
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PSD by parameterizing it as a power law, P(f) oc f~ a , and 
quantifying it by the slope index a which was found to be 
between 1 and 2 for most of Seyfert 1 galaxies (e.g. Lawrence 
& Papadakis 1993; Edelson & Nandra 1999; Markowitz et 
al. 2003). 

A weakness of the Fourier PSD approach to highly and 
aperiodically variable time-series of AGN lies in its limited 
information that can be recovered. Since we are working only 
with a signal decomposition into frequency components, it is 
difficult to dig up an idea about e.g. their life-times and fre- 
quency distributions that are responsible for the observed 
PSD shape. Vaughan et al. (2003) discussed fundamental 
tools for examination of AGN light curve properties along 
the time, however, they skipped the discussion on the meth- 
ods pertaining to the studies of time-frequency domain. So 
far, minor efforts have been undertaken in that direction 
for X-ray accretion astrophysics. A promising tool to scan 
the local and shortly lasting quasi-periodic fluctuations is 
associated with an application of the wavelet analysis (see 
e.g. Lachowicz & Czerny 2005, hereafter LC05, and refer- 
ences therein) . It allows for a signal decomposition into two- 
dimensional space of parameters and to gather a knowledge 
about their deployment. 

Contrary to the results obtained for Galactic black-hole 
binary system of Cyg X-l (LC05), the usage of wavelet trans- 
form to X-ray data of AGN constitutes a challenging step 
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Figure 1. 0.5-10 keV EPIC-pn light curves of MCG-6-30-15 
analyzed in this paper with the wavelet analysis. Here, for visual 
purposes, data were plotted with 100 s resolution. 
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Figure 2. 0.5-10 keV EPIC-pn light curves of NGC 4051 used 
in the wavelet analysis (plotted with 100 s resolution). 



analysis issue stressing some crucial properties. The results 
of data analysis and their discussion are presented in Section 
4. We conclude in Section 5. 



towards understanding their time variability. As a subject 
of our study we choose MCG-6-30-15 and NGC 4051 which 
are frequently observed sources with a high count rate in 
X-ray energy band. 

MCG-6-30-15 is a nearby Narrow-Line Seyfert 1 galaxy 
located at z = 0.007749. It belongs to the one of the bright- 
est in its own class. It shows up a very broad FeK line 
what might suggest that a central black-hole of a mass 
2.9+i;g x 10 6 Mq (McHardy et al. 2005) is spinning rapidly 
and relativistic effects of space-time affect observed X-ray 
flux properties. A broad study of X-ray variability aspects 
for MCG-6-30-15 have been done by Ponti et al. (2004) [rev. 
108] and Vaughan, Fabian & Nandra (2003) [rev. 301, 302, 
303]. The latter found a typical PSD slope value of a ~ 2.5 
breaking to the flatter one (a w 1) below ~2x 10~ 4 Hz. 
High-resolution PSD revealed a dominant white noise com- 
ponent at frequencies >5 x 10 -3 Hz in 0.2-10 keV band. 
NGC 4051 is another nearby Narrow-Line Seyfert 1 galaxy 
(z — 0.0023) which displays a complex spectral behavior 
and highly variable X-ray flux (e.g. Ponti et al. 2006). Its 
PSD reveals energy independent bending at 8^3 x 10 -4 Hz 
frequency changing a slope from -1.1 to about -2 (McHardy 
et al. 2004). Vaughan & Uttley (2005) showed some evi- 
dences of QPO structure to be present at about 5 x 10 -3 Hz 
which may be either real or arise due to inappropriate con- 
tinuum model used in PSD fitting. A black-hole harbored 
in the heart of NGC 4051 are considered to have a mass of 
2-^11 x 10 5 M Q (Shemmer et al. 2003; McHardy et al. 2004). 

In this paper we apply the wavelet analysis to publicly 
available EPIC-pn/XMM-iVewfon observations of MCG-6- 
30-15 and NGC 4051 in order to examine the time-frequency 
domain of their variability. Our paper is organized as fol- 
lows. In Section 2 we describe the details of data selection 
and reduction. Section 3 contains a discussion of the wavelet 



2 DATA SELECTION AND REDUCTION 

For almost all XMM-Newton (Jansen et al. 2001) observa- 
tions of MCG-6-30-15 and NGC 4051, a pn (Strueder et al. 
2001) and two MOS (Turner et al. 2001) CCD detectors of 
the European Photon Imaging Camera (EPIC) were oper- 
ating simultaneously. Due to increased soft-p + flaring level 
detected in the data coming from the MOS cameras, we de- 
cided to ignore them in our analysis. We processed all odf 
files with the XMM-Newton Science Analysis Software (SAS 
v.7.0) using the most updated calibration files. Pn event lists 
were extracted using the standard processing SAS tool of 
epchain. The source's counts were read out from a circular 
region (centered for MCG-6-30-15 at R.A. = 13 h 35 m 53=8, 
8 = -34° 17' 44" and R.A. = 12 h 3 m 9 ! :6, 8 = +44°31'53" for 
NGC 4051) of the radius equal 36" (rev.301) and 40" for 
the rest. A SAS tool tabgtigen was applied to create good 
time interval (GTI) masks which were used in filtering the 
original pn event files. The background scientific products 
were extracted from the circular regions (of the same ra- 
dius as used for the source) on the same chip-set around the 
galaxies. The analysis of the pn data was performed using 
both single and double events (i.e. pattern < 4). In result, 
0.5-10 keV light curve with time resolution of AT = 10 s for 
each XMM-Newton data set was extracted and corrected for 
the background (Figureri]and|2j- We choose data analysis in 
that particular energy band to have highest possible signal- 
to-noise ratio 1 . In addition, to work only with evenly dis- 
tributed records, we applied an interpolation of that counts 
where X-ray flux: (a) dropped to zero, and/or (b) lay above 
3a around the mean (within a given segment; see below) 
and (c) a number of such nearby laid points was less then 

1 Energy-dependent wavelet analysis is the subject of our upcom- 
ing paper (Lachowicz et al., in preparation). 
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10. That correction comprised not more than 1 per cent of 
all data points in all seven data sets. Periods of distorted 
flux level (i.e. close to zero), continuous and longer than 200 
s were skipped. 

For the purpose of the wavelet analysis, each resulting 



light curve was divided into T sf 



10 s long segments 



what establishes a sample of data to be analyzed. Their total 
number, M scg , yielded 41 and 13 for MCG-6-30-15 and NGC 
4051, respectively. 

All performed calculations of the continuous wavelet 
transform were done using the matlab software provided 
by Torrence & Compo (1998, hereafter TC98) and Grinsted, 
Moore & Jevrejeva (2004). 



3 WAVELET ANALYSIS 

The wavelet analysis constitutes an alternative tool in a 
time-series analysis when confronted with the classical ap- 
proach based on the Fourier techniques. In the latter, a sig- 
nal decomposition with the analyzing function, ip, of the 
form of e 1 -' ' allows to detect strictly periodic modulations 
buried in the signal. Signal integration in the Fourier trans- 
form is performed along its total duration. Thus, one cannot 
claim anything about time-frequency properties of periodic 
components like their life-times or time-localizations. A sim- 
ple solution to this problem has been found via the appli- 
cation of windowed Fourier transform (e.g. Cohen 1995), 
i.e. computation performed for a sliding segment of time- 
series (much shorter than a whole signal span) along the 
time axis. That returns a time-frequency chessboard of equal 
resolution everywhere dictated by the Heisenberg-Gabor un- 
certainty principle (Weyl 1928; Gabor 1946; Brillouin 1959; 
Flandrin 1999). A drawback of that method lies in the loss 
of information about fine and temporal quasi-periodic fluc- 
tuations (for example at high frequencies) which are simply 
smeared due to insufficient resolution. 

The wavelet analysis treats this problem more friendly 
as a basic concept standing for its usefulness relay on probing 
time- frequency plane at different frequencies with different 
resolutions. It can be achieved by selecting the analyzing 
function tp of finite duration in time that meets, in gen- 
eral, J 00 ip{t)dt — condition. Making a scaling and time- 
shifting of ip, at selected time duration of the light curve, 
one can reach and examine a desired field of time-frequency 
space (see below). 

Talking about wavelet analysis we mainly refer to the 
computation of a wavelet power spectrum defined as the 
normalized square of the modulus of the wavelet transform: 



W = £\w k (a m )\ 



(1) 



Otherwise we stress the term we mean, for example, cross- 
wavelet spectrum or wavelet coherence (Kelly et al. 2003; 
Grinsted, Moore & Jevrejeva 2004). 

Here, for an introduction to the concept of continuous 
wavelet transform, we refer a reader to our first paper of 
LC05 and references therein. We find there the best ana- 
lyzing function (a mother wavelet) for X-ray signals to be 

2 Hereafter, we refer to each MCG-6-30-15 and NGC 4051 light 
curve by its XMM-Newton revolution number at which the galaxy 
was observed, e.g. rev. 301. 



Figure 3. A chessboard of time-frequency signal decomposition 
in case of wavelet transform. Grey-color rilled areas represent the 
exemplary non-zero coefficients of the wavelet decomposition with 
linearly increasing frequencies. 



the Morlet wavelet, ip(t) = 7 r" 1/4 e i2 ' r/ot e~ t2/2 , with 2tt/ 
parameter set to 6. It oscillates due to a term of e 1 *, has a 
complex form and suits perfectly to probe a quasi-periodic 
modulations at different frequencies. 

Below, we extend our previous discussion on wavelet 
analysis (LC05) paying a special attention to its usage in 
case of astronomical time-series following the Poisson distri- 
bution and to general wavelet map characteristics. 



3.1 Heisenberg-Gabor uncertainty principle 

Despite a promising step forward, the wavelet analysis suf- 
fers from the same problem encountered in the Fourier tech- 
niques, namely, the Heisenberg-Gabor uncertainty principle. 
Applied to time-series it says that one cannot determine in 
an exact way both time and frequency location of extracted 
features (wavelet peaks) simultaneously. Certain uncertainty 
a exists and the product of frequency band-width Af with a 
time duration of the signal At cannot be less than a minimal 
value given by: 



g = AiA/ > — 

~~ 47T 



(2) 



(e.g. Flandrin 1999). One can represent the spread of the 
wavelet in the time-frequency domain by drawing boxes of 
side lengths 2Ai by 2A/ known as Heisenberg boxes. For the 
Gaussian windowed functions as used in the Morlet wavelet, 
the product of AtAj is exactly equal 1/4-7T. 

An exemplary chessboard of time-frequency signal de- 
composition for the wavelet transform is shown in Figure [2] 
Shapes of the boxes represent the actual resolution at the 
wavelet map in a function of linearly increasing frequencies. 
At high frequencies an increase in spread of spectral com- 
ponents is related to decrease in a wavelet width in a time 
domain. A reverse situation takes place at lower frequencies. 
Due to this fact, a sinusoidal fluctuation of the finite and the 
same time duration, located at frequencies fi,fo,fa where 
fi < H < fs, will differ in their time- frequency properties 
in the wavelet map (Figure [3J. We will return to this point 
in Section f3. 41 
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3.2 Significance of wavelet peaks for Poisson 
signals 

A map of wavelet power spectrum shows up a distribution 
of spectral components in the time-frequency plane. Single 
and finite oscillation, of sufficiently high amplitude to be 
detected, marks itself in the map as a (wavelet) peak. To 
every peak a certain statistical significance can be assessed. 
It was done for the first time by TC98. We say that a wavelet 
peak is significant with an assumed per cent confidence when 
it is above a certain background spectrum. The latter is 
determined by the mean power spectrum of analyzed time- 
series. 

TC98 considered and adopted in the publicly available 
software the background power spectrum of the following 
form: 



(3) 



II - ae-^'l 2 



which refers to an assumption that tested signal can be con- 
sidered as a realization of univariate lag-1 autoregressive 
(AR1) process: 



x\ 



axi-i 



(4) 



where a is lag-1 autocorrelation, zi denotes a random vari- 
able drawn from the Gaussian distribution and xq = 0. 
They showed that the distribution of the local wavelet power 
spectrum follows xi distribution with v degrees of freedom 
(equal 1 or 2 for a real or complex wavelet, respectively). 

However, a simple Lorentzian-shape power spectrum 
given by equation © is not a good representation of a 
broad-band PSD shape of MCG-6-30-15 and NGC 4051 (e.g. 
Vaughan & Uttley 2005); the latters are better fitted with 
broken power-law models plus white noise. 

In the present paper we aim at modeling very short light 
curves (T seg ~ 10 4 s, see Section |2J which are dominated by 
the Poisson white noise component. Therefore, in order to 
define the background spectrum for a signal of Poisson dis- 
tribution, we recall equation A7 of LC05 that defines, here, a 
normalized square of modulus of discrete form of continuous 
wavelet transform in terms of the Fourier transform as: 



£\wk(a m )\ 2 = 



At 



1/2 ' 



2=1 



#(27! 



(5) 



where the discrete Fourier transform of signal xi is given by 



-i2TTjl/N abs 



(6) 



j denotes a frequency index, JV b s a number of data points 
and x a mean value of Xi 3 . Vaughan et al. (2003) pre- 
sented a discussion on popular normalizations of X-ray 
Fourier power spectra, namely ^4|a;j| 2 . Among them we find 
A — 2Atx~ as the most interesting for our purposes: 

(i) it has the property that for such normalized power spec- 
trum the expected Poisson noise level equals simply 2, and 

(ii) for any signal following the Poisson distribution, result- 
ing power spectrum should be distributed exactly as xli i- e - 
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Figure 4. A simplified wavelet map for simulated 10 sinusoidal 
modulations of the finite time duration (8t = 1500 s) localized 
at frequencies /j = O.OOli Hz (i = 1,2,..., 10). Please note that 
frequency axis is plotted logarithmically contrary to Figure|!J] Sec 
Section 13.41 for details. 



with 2 degrees of freedom (Leahy et al. 1983). That fact 
allows us thus to construct the background spectrum for 
Poisson data as follows: 



pPoiM = 2 



(7) 

assuming that the input light curve xi is in units of cts s _1 
and, what most important, the normalization factor of the 
wavelet power spectrum equals: 

£ = AN obs = 2Atx~ 1 . (8) 

Therefore the integrated wavelet spectrum over time (global 
wavelet spectrum) defined as previously (LC05) in the form 
of: 



G{a m ) = 



k 



\wk{a m )\ 



(9) 



the same as A | a; j | 2 . 



will hold the units of cts s Hz~ 
Assuming normalization of JSJ and recalling property (ii) 
given above, by the analogy to TC98, the local wavelet power 
spectrum ought to be \\ distributed, i.e. 



£,\wk(a m )\ 



1 nPoiss 2 

2^ X2 



(10) 



for each time location k and at particular scale a m . In gen- 
eral, this is always provided by the central limit theorem: the 
Fourier transform is the sum over time of a function of all xi 
and in the limit of large iV bs this sum tends asymptotically 
to Gaussian distribution. 

Summarizing the above, for shortly spanned X-ray light 
curve (provided widely e.g. by the detectors on board XMM- 
Newton or RXTE facilities) any wavelet peak standing sig- 
nificantly above the mean background spectrum given by Q 
can be considered as a true feature at assumed significance 
level 4 . 



We note that a factor N ohs in equation A6 of LC05 appeared 
by mistake. 



4 Fourier spectrum of the Poisson process is frequency indepen- 
dent, thus pP° lss does not depend on j. Here, we neglect a red 
noise component since we focus at high frequencies. 
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Figure 5. Upper: a dependence of frequency spread of the 
wavelet peak in a function of its location at frequency / in the 
wavelet map. Bottom: a dependence of ratio of observed peak du- 
ration to the actual one in a function of signal localization at a 
given frequency (valid for AT = 10 s; see Section l3.4l . 

3.3 Statistical description of a wavelet map 

Depending on the frequency band that we wish to analyze 
with the wavelet analysis, a long-lasting light curve (of du- 
ration T and time resolution AT) is advised to be divided 
into M sog segments (of duration T sog ). Therefore, for ev- 
ery segment the wavelet power spectrum can be computed 
individually. In the resulting wavelet maps, as the wavelet 
peak one may define the contour of its statistical significance 
calculated at the given level. In order to describe quantita- 
tively the content of peaks and their time-frequency proper- 
ties along T, we propose: 

(i) to determine a median position of each peak: in time, 
[to] , and in frequency, |_/ol ; 

(ii) to give the total number of peaks, N p , per light curve; 

(iii) to determine the average number of peaks per seg- 
ment, N p /M seg ; 

(iv) to read out a horizontal and vertical span and denote 
them as observed peak duration <5£ b s and its width in fre- 
quency <5/obs, respectively, where 8t b s should be corrected 
for the effect of the peak location at frequency / (see Section 
EH for details); 

(v) to calculate a median value of (sinusoidal) oscillations, 
L-Wosc] , i.e. for case when assumed ip(t) is the Morlet wavelet 
and thus N osc = L/o~| <5i for each peak. 

The above set of parameters can be treated as a basic one 
and used for the future purposes of comparison of results 
from the wavelet analysis between different X-ray sources. 

3.4 Wavelet peak resolution as a function of 
frequency 

In order to quantify the dependencies of the signal duration 
St and its width in frequency domain, <5/, as a function of the 
peak location at frequency /, the following simulation was 
performed. We generated a light curve of the total duration 
of T = 40000 s at bin time AT = 10 s with 10 sine function 
yi = sin(2n fit) + const, fi — O.OOli and i = 1,2, 10. Each 



%ji was assumed to last through a constant period of time of 
8t — 1500 s, separated from each other by about 1500-2000 
s in such a way that const value was assumed between the 
sinusoidal parts of the signal. No noise component was added 
as we were interested only in a basic response of the wavelet 
transform to tested signals. A simplified wavelet map for 
that simulation is presented in Figure [I] 

For a given wavelet peak its time-frequency location in 
the map was determined as given in Section [!0] and marked 
by a dot in the figure. As tested before (LC05), due to irregu- 
larities in the peaks' shapes, in order to obtain best peak po- 
sition, a median value was found to work better rather than a 
mean. Here, it returns a location in frequency with the aver- 
age error of the order of 10~ 4 , namely fi = 0.001 (±0.0001)i 
for i = 2,3, 10. 

In Figure |S] (upper panel) the relation between <5/ b s of 
the peak as a function of / is given. A solid line represents 
the best fit to the data of the quadratic form: 

5 fobs = -18.11/ 2 + 0.39/ + 0.0006. (11) 

Therefore, a quantity of <5/ b s /2 can serve as a rough esti- 
mation of an error we make in attempt of determination of 
the peak position at given frequency. 

A second and non-negligible wavelet map property per- 
tains to the dependence of the ratio of observed peak life- 
time to the actual one (i.e. assumed or/and known in ad- 
vance) as a function of /. It is not constant as one might 
expect. Figure^] (bottom panel) presents this relation. Here, 
a solid line denotes the best fit of assumed model y(f) = 
a/f + b. Therefore, one can correct the observed life-time 
of the wavelet peak, Stobs, for its real duration, 8t, in the 
following way: 

5t = 5t ohs (^j^- + Q.77\ . (12) 

Formula 111 21 can be treated as a good approximation of St 
with a relative error about 2 per cent for the wavelet map 
calculation in 0.0008-0.01 Hz band. 

We would like also to note that relations li lt and 1121 
ought to be recalculated depending on our choice of AT. 
Here, they are valid for AT = 10 s. 



4 RESULTS AND DISCUSSION 

We construct the wavelet power spectrum maps for M seg 
segments of data and characterize the source activity dur- 
ing a given period by extracted wavelet peaks at the signifi- 
cance level of 98 per cent 5 . Figure |S| (upper panel) presents 
an exemplary outlook of such time-frequency decomposition 
performed for one segment of data (rev. 301). In the bottom 
panel of the same figure we show the corresponding Fourier 
power spectrum calculated for the same frequency band as 
given in the wavelet map. The red straight line represents 
the mean background spectrum expected to be the Poisson 
noise level (Section 13. 21 . For every light curve we calculate 
its statistical characteristics as described in Section f3.3l and 



5 We explain our choice of that particular significance level in 
Section IO 
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Table 1. Properties of the localized peaks in the wavelet maps between 0.002-0.007 Hz (0.5-10 keV) for the assumed significance level 
of 98 per cent. M Beg denotes the number of segments the given light curve was divided into. N p stands for the total number of detected 
wavelet peaks. [St] and [<5/l are the median values of wavelet peak life-times and frequency widths, respectively, whereas [Aosc] is the 
median value of the sinusoidal oscillations among all peaks, x denotes the mean count rate of the light curve. 



Rev. M scg N p N p /M scg [St] [Sf\ [N oac ] x 

(s) (Hz) (cts s" 1 ) 
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Figure 6. Upper panel: exemplary wavelet map calculated for 
10 4 s long segment of MCG-6-30-15 light curve (rev. 301). Top 
values of wavelet power are denoted by gradual brightening of 
the color. White solid contours denote significance level of 98 per 
cent for extracted peaks. A cone of influence (COI) is marked as 
a white area at both sides of the map and corresponds to that 
space of wavelet spectrum for which the results are affected by 
the effect of signal finite duration. Bottom: corresponding Fourier 
power spectrum. Red straight line denotes the mean background 
spectrum at the expected Poisson noise level. 

summarize them in Table In Table |2] we compare glob- 
ally the properties between MCG-6-30-15, NGC 4051 and 
simulated light curves of Poisson noise. 

4.1 Statistical properties of wavelet activity 

Table [5] reveals that the median life-times of peaks for each 
object remain very similar, however, they change within 
short period of time (Table 0. In particular, for MCG-6- 
30-15 observed 3 times during consecutive 5 days (rev. 301- 
303), the maximal difference in [St] yields about 213 s. More 



dramatic changes we can see for NGC 4051, i.e. of about 
350 s, but, in this case, separated by 18 months in the source 
observation moments. That may indicate that character of 
physical processes operating in the source at time-scales be- 
tween 142 s and 500 s is not constant and may change even 
by a factor of ~2 or more for [St] within a few days. 

Interestingly, the average life-time of the peaks does not 
seem to be constant function of frequency at which they 
occur. It is shown in the upper panel of Figure[7| For MCG- 
6-30-15 the existing variability structures tend to live longer 
at longer time-scales but for NGC 4051 the same behavior is 
different. Oscillations live longer with increasing frequency, 
peaking at about 4.5 x 10~ 3 Hz and then suddenly drop down 
by a factor of 20. The cause of this rapid drop is unclear at 
this stage of research. 

The same cannot be read out for the ratio of the peak 
frequency to its width, denoted widely in the literature as 
a quality factor, Q — f/Sf. For both sources is remains 
approximately constant with the average value of (Q) ~ 3. 

The results obtained from the wavelet analysis also pro- 
vides us with the opportunity of the histogram representa- 
tion (hereafter also referred to as the wavelet power his- 
togram) of the total number of detected oscillations at a 
given frequency. It can be considered as an analogue to the 
global wavelet power spectrum or the Fourier power spec- 
trum. We performed a computation of such histograms for 
every light curve separately and displayed the result in the 
form of stacked bar graph for MCG-6-30-15 (Figure|SJi) and 
NGC 4051 (Figure [Hp), respectively. These graphs show up 
contributing amounts of the wavelet peaks detected in each 
observation to their overall activity at given frequency. As 
denoted by Ponti et al. (2006), the second XMM-Newton 
observation of NGC 4051 (rev. 541) was timed to coincide 
with an extended period of its low X-ray emission that is 
why the total number of extracted events from this obser- 
vation is very low (Table 0. 

For MCG-6-30-15 the histogram reveals a sort of bi- 
modal distribution with a distinctive peak at ~2.5 x 10 -3 Hz 
and a broad one centered at ~5 x 10~ 3 Hz, respectively. It 
is difficult to judge its reality looking only at the Fourier 
power spectra (cf. figure 4 in Vaughan, Fabian & Nandra 
2003). However, it is worth to mark that our high frequency 
(HF) peak stays in a good agreement with some predictions 
for HF QPO (to be in 0.5 - 3 x 10~ 3 Hz) of Vaughan & 
Uttley (2005). A less pronounced bimodal character can be 



© 2006 RAS, MNRAS 000, 000-000 



Wavelet analysis of MCG-6-30-15 and NGC 4051 7 



Table 2. Comparison of the wavelet peak characteristics ob- 
tained from the analysis of MCG-6-30-15 and NGC 4051 data 
and simulated time-series containing only a pure Poisson noise. 
Rows (up to down) provide information on: number of segments, 
total number of extracted wavelet peaks, the average number of 
peaks per segment, the median life-time of peaks, the median 
width of peaks and the median number of sinusoidal oscillations 
within duration of the peak, respectively. See Section l4.2l for de- 
tails. 



Quantity 


MCG-6-30-15 


NGC 4051 


Simulation 




41 


13 


1000 




129 


84 


1553 


N p /M seg 


3.14 


6.46 


1.55 


[St] 


366.6 s 


385.6 s 


308.5 s 


W 


0.0012 Hz 


0.0012 Hz 


0.0011 Hz 


[N OBC ] 


1.56 


1.48 


1.33 




O - 




Frequency, f [Hz] 



Figure 7. A relation between (upper panel) the average life- 
time of wavelet peaks, (St), and (bottom part) the average quality 
factor, (Q), of detected wavelet peaks as a function of frequency 
for MCG-6-30-15 (open circles) and NGC 4051 (filled diamonds). 




«N p » = 6.17 


(c) 

































^^^^^^^^^^^^^^^^^^ m 
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Frequency, f [Hz] x1CT 3 

Figure 8. Histogram of wavelet peak distribution in a function 
of their location in frequency: (a) for MCG-6-30-15 where colors 
refer to the single histogram values for each light curve separately 
(see legend) whereas joined together represent a contribution to 
the total peak activity in a given frequency band, (b) for NGC 
4051, and (c) for simulated time-series of the pure Poisson noise 
(see Section l4.2l for details). 



seen for the NGC 4051 wavelet power histogram. Here, an 
indication of QPO at 5 x 10~ 3 Hz in NGC 4051 (mentioned 
in Section differs comparing to our results. 

Histograms for both galaxies were calculated at as- 
sumed significance level of 98 per cent. We have checked the 
results for 95-99 per cent band of significance and we have 
found that only for 99 per cent level the histograms failed 
to show up a bimodal distribution due to lowest number of 
extracted peaks. Therefore, 98 per cent level was assumed 
by us as the highest possible level in order to perform cal- 
culations and the analysis of results. 

4.2 Wavelet analysis of pure Poisson noise and an 
excess wavelet power histogram 

Vaughan, Fabian & Nandra (2003) performed a detailed 
Fourier power spectrum analysis for MCG-6-30-15 obser- 
vations at rev. 301, 302 and 303. For high-frequency pe- 
riodograms estimated using AT = 10 s they found that 
a red-noise PSD of the source can be seen for frequencies 
<3x 10~ 3 Hz. It is intuitively confirmed in the bottom part 



of Figure [S] though, in this case T sog ~ T/8, i.e. time-series 
is much shorter than a whole light curve duration. There- 
fore, based on the Fourier technique results one can con- 
sider 0.003-0.007 Hz fully dominated by the Poisson noise. 
Hereby, we challenge this statement to may be not fully cor- 
rect by performing the following Monte-Carlo simulation. 

To find out about the ratio of all wavelet peaks ap- 
pearing in the wavelet maps due to pure Poisson noise to 
the total number of peaks observed for MCG-6-30-15 (and 
NGC 4051), we generated 1000 light curves of duration 10 4 s 
each at AT = 10 s. Here, we followed the same procedure 
as described in section 5.1 of LC05 (equation 13), however 
assuming a constant value equal to the average count rate in 
MCG-6-30-15 and NGC 4051 EPIC-pn data, i.e. ~15 cts s" 1 
(0.5-10 keV). For every single time-series the wavelet power 
spectrum in 0.002-0.007 Hz frequency band was calculated 
and statistical characteristics (of Section f3,3H were done. 

Table confronts the wavelet peak properties between 
these obtained for two galaxies (all light curves included) 
and from the simulation. The most striking result pertains 
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4 5 
Frequency, f [Hz] 



x 10"' 



rev. 1 08 
+301 
+302 
+303 



rev. 263 
+541 



Figure 9. Excess wavelet power histogram for combined (a) 
MCG-6-30-15 and (b) NGC 4051 light curves, respectively. See 
Section l4.2l for details. 




x 10 



Figure 10. Frequency-frequency plots for twin peak QPOs in 
MCG-6-30-15 and NGC 4051. The plots are remarkably similar. 
They suggest that two resonances may be present here, the 3:2 
resonance that is also common in the Galactic black hole and 
neutron star sources, and a 2:1 resonance. 



to the ratio of the average number of peaks (per segment) 
extracted for MCG-6-30-15 and for pure Poisson process 
realization 6 which is equal to ~2. Therefore, this ratio says 
that statistically about half of all observed wavelet peaks 
in 0.5-10 keV X-ray MCG-6-30-15 emission is of the pure 
noise origin. NGC 4051 occurs to be more active with the 
average number of peaks per segment of ~6.5 where only 24 
per cent of them can arise due to noise. 

In order to visualize the distribution of wavelet peaks 
being a pure noise contamination, we construct the noise 
histograms in the following way. We group m x M seg sim- 
ulated light curves of Poisson noise into m parts. For each 
part (composed of M scg segments) first we sort the list of 
peaks according to frequency and then we calculate a his- 
togram (10 bins with 0.0005 Hz bin-width). Finally, for all 
of m histograms, in a given histogram-bin we compute the 
mean number of peaks, {N p ), and associated error. The lat- 
ter appears to be of the order of 10~ 2 . Here, m was assumed 
to be equal 24 and 76 for MCG-6-30-15 and NGC 4051, re- 
spectively. In result we found that both histograms revealed 
a uniform distribution of peaks as a function of frequency 
with the mean value, {(N p )) = 6.17 and 1.96 for m — 24 and 
76, respectively. Figure^ shows up the noise histogram cor- 
responding to m = 24 and A/ scg = 41 (the error bars have 
been omitted). 

This outcome can serve us as a good estimation of the 
average number of wavelet peaks in the light-curves that 
are of a pure noise origin. Therefore, it also allows us to 
subtract ({N p }) value from the data (Figure |Hp,b) obtaining 
an excess wavelet power histogram presented in Figure |U] for 
both galaxies. Thus, the excess histogram provides us with 
information on the wavelet peak distribution in the source 
disentangled from the average level of wavelet peaks arising 
only due to noise process. 

6 All detected wavelet peaks for Poisson noise signals correspond 
to real and significant wavelet powers standing out above the 
mean background spectrum. 



4.3 The frequency-frequency (/high vs. /i ow ) plot 

Quasi-bimodal distribution of peaks that can be seen in 
Figure [5^ suggests that MCG-6-30-15 may show double 
peak QPOs (with frequencies /high and /i ow ) as indeed 
many Galactic black hole and neutron sources do. One 
can read out from this figure that the ratio is equal about 
/high : /low ^2:1 assuming /high and /i ow to be ~0.0055 Hz 
and ~0.0025 Hz, respectively. If the histogram peak at 
~0. 00375 Hz goes into a play that indicates also at 3:2 fre- 
quency ratio. Similar ratios can be found for the excess his- 
togram of NGC 4051 (Figure Ep). 

We decided to investigate this case more carefully. We 
divided the list of wavelet peaks detected in each of four 
MCG-6-30-15 light curves into two parts of equal amounts. 
Thus, for 8 parts of data we calculated 8 histograms (10 bins) 
of the peak distribution as a function of frequency. Each 
histogram revealed two distinctive and well separated peaks 
at frequencies denoted here as /i ow and /high, respectively. In 
order to increase the precision of frequency determination we 
also checked the histograms with 11 and 12 bin resolution. 
Where needed, 8 and 9 bin widths were also investigated. 
By the average error of determination of peak position we 
assumed a half of the peak width. In result, we plotted the 
relation of /high as a function of /i ow for given 8 pairs in one 
figure (see Figure ITTJ1 . We have repeated the same procedure 
for NGC 4051 (5 pairs) obtaining very similar results. 

Rather remarkably, for both sources a few points locate 
along a 3:2 frequency ratio line, and the rest close to a 2:1 
line. No other combinations are observed. 

4.4 Double peak QPOs in MCG-6-30-15 and 
NGC 4051: a link to Galactic kHz QPOs? 

Results discussed in the previous Section and displayed in 
Figure fTTJl suggest that double peak QPOs in both MCG-6- 
30-15 and NGC 4051 have commensurable frequency ratios, 
3:2 and 2:1. Abramowicz & Kluzniak (2001) first noticed 
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Figure 11. McClintock & Rcmillard (2003) scaling for micro- 
quasars extended down to AGN region with marked results found 
for MCG-6-30-15 and NGC 4051 within our wavelet research. If 
real, it may be considered as the link between QPOs observed in 
AGNs and kHz QPOs seen in X-ray binary systems (see Section 
I4,4l for discussion). 



the 3:2 ratio (and stressed its importance) in the X-ray vari- 
ability data of the microquasar GRO J1655-40 analyzed by 
Strohmayer (2001); see also McClintock & Remillard (2003). 
Abramowicz et al. (2003) argued that the 3:2 ratio is also 
evident in the case of the neutron star QPOs. Although the 
ratio is varying, it clearly distinguishes the 3:2 value (see Bel- 
loni et al. 2005 for a criticism). There is a consensus today 
that the commensurable frequency ratios, and 3:2 in partic- 
ular, are an important property of the QPO phenomenon in 
microquasars and neutron star sources (for recent references 
see Torok et al. 2005). It has been pointed out by Kluzniak 
& Abramowicz (2000) that a commensurable frequency ratio 
suggests a resonance, and that the double peak QPO in mi- 
croquasars and neutron stars are due to a similar non-linear 
resonance between two modes of accretion disk oscillations 
in strong gravity. 

Our results may suggest that the same 3:2 resonance 
operates also in the MCG-6-30-15 and NGC 4051 accre- 
tion discs. The full theoretical discussion of this suggestion 
will be presented in a follow-up paper, now in preparation. 
Here, we would like to make only one comment, on the QPO 
frequency scaling with the inverse mass of the source, that 
is theoretically expected in strong gravity. McClintock & 
Remillard (2001) noticed that for three microquasars with 
known mass (see also Torok et al. 2005) , the QPO frequency 
indeed scales according to: 



(13) 



In Figure [TTI we show that this scaling holds for NGC 
4051 with the commonly adopted mass M = 2-f- 11 x 10 s M 
(Shemmer et al. 2003; McHardy et al. 2004), but for MCG- 
6-30-15 with a rather low mass estimate, M ~ 5 x 10 5 Mq 
that according to McHardy et al. (2005) is possible, but not 
very likely. 



5 CONCLUSIONS 

We have presented the results of the wavelet analysis for two 
Narrow-Line Syl galaxies, MCG-6-30-15 and NGC 4051, in 
0.5-10 keV energy band. The main points of our research 
can be summarized as follows: 

(i) We explored the observational properties of time- 
frequency wavelet chessboard (Section |3J and proposed a 
set of common parameters that can be used in the future in 
order to compare wavelet power spectrum results obtained 
for different X-ray sources (Section l3.31 : 

(ii) We discussed and showed (Section 13. 21 that the most 
proper mean background spectrum required in the tests for 
the significance of wavelet peaks in case of Poisson signals 
can be assumed as given by equation Q; 

(iii) We showed that the wavelet (histogram) analysis is a 
powerful and efficient method to extract real signals buried 
in the Poisson noise which are present in the data for very 
limited fraction of time, and how to disentangle them from 
that noise (Section 14. 2H : 

(iv) Finally, we seem to detect pairs of QPO oscillations 
in MCG-6-30-15 and NGC 4051, which are either in 2:1 or 
3:2 resonance (Section 14.31 . If confirmed in the future with 
extended number of observations for the same and/or other 
bright AGNs, our finding may indicate on a link between 
QPOs observed in AGNs and kHz QPOs seen in X-ray bi- 
nary systems (Section I4.4|l . 
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